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Inverse Compton X-rays from Giant Radio Galaxies at z ~1 
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ABSTRACT 

We report XMM-Newton observations of three FR II radio galaxies at redshifts between 0.85 
and 1.34, which show extended diffuse X-ray emission within the radio lobes, likely due to 
inverse-Compton up-scattering of the cosmic microwave background. Under this assumption, 
through spectrum- fitting together with archival VLA radio observations, we derive an inde- 
pendent estimate of the magnetic field in the radio lobes of 3C 469.1 and compare it with the 
equipartition value. We find concordance between these two estimates as long as the turnover 
in the energy distribution of the particles occurs at a Lorentz factor in excess of ~ 250. We de- 
termine the total energy in relativistic particles in the radio emitting lobes of all three sources 
to range between 3 — 8 x 10^^ erg. The nuclei of these X-ray sources are heavily-absorbed 
powerful AGN. 

Key words: X-rays: galaxies — galaxies: jets — galaxies: high-redshift — X-rays: individual 
(3C 469.1, MRC 2216-206, MRC 0947-249). 



1 INTRODUCTION 

Emission from inverse-Compton scattered CMB (ICCMB) pho- 
tons, in the form of diffuse, extended X-rays between the nucleus 
and radio hot spots, has been det ected in a number of radio galaxies 
at c osmological distances (e.g. 'Fabian et al. 20031 : iBlundell et all 
[2OO6 : Erlund et al. 2008; Johnson et al. 2001) and at low redshift 
(e.g. ICroston et all 1200 5). ICCMB is a tracer of old, spent syn- 
chrotron plasma and its usefulness as a diagnostic of the magnetic 
field and electron energy distributio n in AGN lobes s tems from the 
fact that it is not redshift-dimmed (ISchwartj|2002h . unlike other 
sources of continuum X-ray emission. 

Together with radio observations, ICCMB X-ray fluxes can be 
used to constrain the local magnetic field e xperienced by the radio- 
emitting plasma in the lobes of giant FR II (iFanaroff & Rilevll974h 
radio sources (see e.g. Erlund et al. 2006). Previous observations 
of ICCMB X-rays from radio galaxies have found deviations of 
the observed m agnetic field from p redictions based on equipartition 
arguments (see lCroston et al.|[2005 , and references therein). While 
this can be r esolved if the energy distribution turns over at the low- 
energy end dBlundell et al.l l2006) at relatively high Lorentz factors 
(7min~ 10^), estimates of the low-energy turn-over of the electron 
distribution in all observed radio galaxies and quasars have been 
largely inconclusive, with no consistent picture amongst the entire 
population of radio sources having yet emerged. 

We present XMM-Newton observations of the FR II radio 
galaxy 3C 469.1, where we have detected extended X-ray emis- 
sion, which is likely ICCMB. We extract and fit the X-ray spec- 



trum for this source and show the photon spectral index, F to lie 
between 0.9 and 2.1, consistent with the radio synchrotron spectral 
index of 1.97 determined from NEE0 archival radio photometry. 
We use the calculated unabsorbed X-ray luminosity along with a 
new look at archival 1.5 GHz VLA data for this source to con- 
strain the magnetic field, without assumptions of equipartition. Fi- 
nally, we present X-ray observations of two more giant radio galax- 
ies, MRC 2216-206 and MRC 0947-249, where ICCMB emission 
seems to be also apparent. 

3C 469.1, MRC 2216-206 and MRC 0947-249 are located at 
redshifts of 1.336, 1.148 and 0.854 respectively, corresponding to 
linear scales of 8.6 kpc/^ 8.4 kpc/^and 7.8 kpc/^ in our adopted 
cosmology of Hq — 70 km s~^Mpc~^ and Qa = 0.73 and as- 
suming a flat universe. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 X-ray and radio imaging 

Our XMM-Newton observations of 3C 469.1 were taken on 2008 
September 29 and consist of 17.4 ks of EPIC-PN and 21.3 ks of 
EPIC-MOS data. Although the MOS data give roughly consistent 
results, they are much noisier and are not presented here. The EPIC- 
PN data were reduced using the standard XMM-SAS pipeline tasks 
EPCHAIN to give 443 ± 21 counts (of which 110 ± 4 are back- 
ground) in the 0.3-10 keV band. 

The source clearly shows extended emission at energies below 
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1.2 keV (Fig.[T]). We estimate the detection significance of the en- 
tire extended emission to be 10-cr and of the northernmost region of 
the quasi-Unear structure to be 5-cr. Above 2 keV, a central core re- 
mains, which is consistent with being a point source. The extended 
linear morphology of the soft X-ray emission precludes it being 
due to hot cluster gas. 

Radio data from all three targets came from the VLA archive. 
In the case of 3C 469.1, we used data from the A, B and C config- 
urations (originally published by Kharb et al. 2008) to sample all 
the structure; for our other two targets only B -configuration data 
were available. Interference was excised from the data and 3C48 
was used to set the absolute flux calibration; a number of rounds 
of phase self-calibration were performed, before a final iteration of 
phase and amplitude self-calibration. The final self-calibrated 1.5 
GHz radio image for 3C 469.1 is shown in Figure[T] The cospatial 
nature of the radio synchrotron emission compared with the 
soft X-rays from this source suggests a fundamental connection 
between the emitting particle populations within the lobes. 

The flux density of the extended radio emission, after subtract- 
ing the contribution from the hotspots, is about 300 mJy, with the 
North and South radio lobes contributing 0.256 Jy and 0.046 Jy, 
respectively. We note that the southern lobe is not detected in the 
soft X-ray image. A possible explanation for the extended emis- 
sion is Inver se-Compton scatteri ng of infrared photons from the 
host galaxy (iBrunetti et al.ll 1997b into the line of sight by plasma 
within the lobes, which would produce more intense X-rays in the 
lobe directed away from the observer. However, in order for the en- 
ergy density of the nuclear photons to be comparable to that of the 
microwave background 200 kpc into the lobes, the nucleus would 
be required to have an infrared luminosity of 10^^ erg s~^ , which is 
implausible. Therefore, it is more likely that ICCMB is responsible 
for the observed X-ray flux. The discrepancy in the flux from the 
two lobes remains unresolved. 



2.2 X-ray spectral analysis of 3C 469.1 

The number of counts was insufficient to perform spectral fitting 
of the lobe and core regions separately (an attempt at fitting only 
the extended X-ray emission towards the North with a power law 
yielded a poorly-constrained photon index ofF = 3.2lj g). There- 
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Figure 2. XMM pn spectrum of 3C 469.1. The model fitted is shown in 
Figure[3]and is described by equation [T] with the best fit parameters as shown 
in Table[I] 

fore spectra of the entire source were extracted, with regions on the 
same chip selected as background. The spectrum was binned to 20 
counts bin"^ using GRPPHA and fitted in XSPEC vl2.5.0 (Fig.E]). 

The model used was a set of power laws in photon energy for 
both the extended and core emission, together absorbed photoelec- 
trically by neutral hydrogen within our Galaxy, with allowance for 
additional soft-photon absorption within the host galaxy itself. Due 
to the poor constraints from the data, we assumed the two power 
law indices to be the same. Mathematically, the model of the X-ray 
flux density as a function of the photon energy, E\ 



M{E) = e 



Here Ki and K2 are the normalisations of the emission at 1 keV 
from the lobes and the core, respectively, Fi and F2 are the corre- 
sponding photon indices (set equal), A^h is the absorbing column 
within our Galaxy along the line of sight to the source (1.38 x 10^^ 
cm~^; Kalberla et al. 2005), A^h is the extra absorbing column 
within the host galaxy, (j{E) is the photo-electric cross-section (ig- 
noring Thomson scattering) and z — 1.336 is the redshift. The best 
fit model (x^ = 13.9 for 17 d.o.f.) is shown in Fig.[3]and the best 
fit values for this model are given in Table [T] 
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Figure 3. Best-fit model for XMM pn data for 3C 469.1. The model used 
(equation [TJ is two absorbed power laws with the same photon index, with 
one power-law component further attenuated by an absorbing column (esti- 
mated at 29.4^^^ g X 10~^^ cm~^ at 90 per cent confidence) within the 
host galaxy. 



Parameter Value 



Nominal error 

(90 per cent confidence) 





1.5 


0.6 




29 


±18 




6.3 


+ 1.5 
-1.7 


K2 


24 


+48 
-20 



Table 1. The constraints on the parameters for the best-fitting model (equa- 
tion [T}. A^^ is the extra absorbing column within the host in 10^^ cm~^. 



Ki and K2 are the normalisations of the emission (in 10 
s~^) at 1 keV from the lobes and the core, respectively. 



• keV" 



The best fit model calls for a high absorption column in the 
nuclear component. A plot of the contours of A^h versus V (Fig. 
131) shows the extra absorption to be statistically significant. After 
correcting for Galactic absorption, the total luminosity in the rest 
frame of the source in the 2-10 keV band is 5.75 x 10^^ erg s~^, of 
which 3.25 x lO^'^ erg s~^ comes from the (absorbed) core compo- 
nent. The remaining luminosity of 2.5 x 10^^ erg s~^, correspond- 
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Figure 4. contours at 68 per cent (innermost, black), 90 per cent (middle, 
red) and 99 per cent (outermost, green) confidence between the absorbing 
column in the host galaxy, A^^ and the power-law photon index, V for the 
best fit model for the 3C 469.1 spectrum. 



X-ray bandwidth of '-^ 8 keV) at 2; = 1.336 from the unabsorbed 
(lobe) component, is comparable to the power in the X-ray emitting 
lobes of the giant radio galaxy, 6C0905+39 (1.5 x 10^^ erg s"^ 
lErlund et aLll2008h . After correcting for both Galactic and host ab- 
sorption, the 2 — 10 keV rest-frame luminosity of the nucleus is 
9.8 X 10^^ ergs-\ 

There are 22 photometric data points from 37 MHz to 14.9 
GHz for 3C 469.1 in the NASA Extragalactic Database. A log-log 
plot of flux density versus frequency yields a straight line with a 
slope of —0.97 and a correlation coefficient of — 0.98. This im- 
plies an electron energy distribution proportional to which 
is considerably steeper than the distribution (^~(2^~^) oc E~'^'^) 
implied by the X-ray spectra. However, our X-ray data only poorly 
constrain the photon spectral index and a value of F = 1.97, as 
expected from the radio photometry, is entirely within the (2cr) un- 
certainty (Fig.|4]). Also, this fit does not take into account the core 
photon index separately, which may have a flatter spectrum than 
the lobes, forcing the overall best fit to lower values of F. Further, 
taking into account the noisy MOS data pushes up the value of F 
from 1.5 to 1.7. In what follows, we shall assume an electron en- 
ergy distribution of A^e(^) oc ^"^-^^ (F = 1.97) for this object. 



3 DETERMINATION OF MAGNETIC FIELD 

For an electron population with a power-law distribution of ener- 
gies, Ne{E) oc the synchrotron photons also follow a power- 
law spectrum with a spectral index, a = (n — l)/2. The ratio of 
the flux density from inverse-Compton emissi on at Uc to th e syn- 
chrotron flux density at Vs is given by (see, e.g. lTuckeilll977h 



2.47 X 10"^^ (5.23 x 10^)"^ 



T 

Ik 



B 



a(n) \l Gauss) 



&(n) 



-(c + l) 



(2) 



where B is the magnetic field strength in the radio-emitting lobes 
and T is the temperature of the CMB at the redshift of the 
source, and the constants a{n) and h(n) have b een tabulated by 
Ginzburg & Syrovatsk ii (1965) and Tucked (Il977h . respectively. 

Taking the flux density of the inverse-Compton emission to 
be the X-ray flux density from 2 — 10 keV (rest frame), the flux 
density of the synchrotron emission to be 300 mJy at 1505 MHz, 
a = 0.97, and hvc ~ 3 keV yields a magnetic field of 12.5 /xG. 
For comparison, taking a — 0.7 in the above formula yields B = 
6.8 /iG, keeping all other values the same. It should be noted that 
this analysis requires the observed X-ray and radio fluxes to be co- 
spatial. If the X-ray emission is more extended than the radio, then 
one must use in the above formulae only that part of the X-ray flux 
which overlays the radio and vice-versa. 

One can also determine the magnetic field by assuming the 
energy density in it to be the same as that of the electron population 
(see Longair 1994, eqn. 18.73). Taking the radio emitting lobes of 
3C 469.1 to be a cylindrical region 640 kpc long and 140 kpc in 
diameter and assuming a minimum Lorentz factor of the electron 
population, 7min = 10^ gives an equipartition magnetic field of 
9.0 /iG. While this value is consistent with that determined from 
the X-ray and radio flux densities, the equipartition field depends 
on the volume of the region {B oc fQj- ^ ^ 0.97), which 

is difficult to determine with the low spatial resolution of XMM 
and given the low number of counts. With the assumed volume 
stated above, we must have 7min > 250 for a < 0.97 in order 
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Figure 5. XMM-pn image MRC 22 1 6-206 from . 3 - 1 . 5 ke V (left panel) Figure 6. XMM pn image of MRC 0947-249 from . 3 - 1 . 5 ke V (left panel) 

and 1.5 - 8.0 keV (from right panel), with 4.86 GHz radio contours over- and 1.5 - 8.0 keV (right panel), with 4.86 GHz radio contours overlaid, 

laid. Each panel is 222 arcsec high and 183 arcsec wide. The X-ray images Each panel is 1 14 arcsec high and 95 arcsec wide, 
have been convolved with a Gaussian with FWHM of 3 pixels. 



for the source to be in equipartition and to simultaneously explain 
the ICCMB X-ray flux. The volume emissivity in the ICCMB X- 
rays, 

= 4.19 X 10-^°T^(z)6(n)Ar f ^ T{z) J 

(3) 

integrated over rest frame energies from 2-10 keV and equated 
to the X-ray luminosity of the lobes, yields the number density 
of relativistic electrons ioh^ N/V — 2.2 x 10~^m~^. Integrat- 
ing the electron energies over a power-law distribution with index, 
n between Lorentz factors 7min and 7max now allows for the de- 
termination of the total energy of the electrons responsible for the 
ICCMB. From the relation 

-1/2 



7 ' 



1/2 



B 



10" 



(4) 



the synchrotron emission at 14.9 GHz jLaing & Peacockll98Qh cor- 
responds to Lorentz factors of 3.2 x 10^. Taking 7max to be this 
value and assuming the electron energy spectrum to extend down 
to 7min = 10^, the total energy in the plasma between 7min and 



1^ 



is 
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1 X 10 erg , which is comparable to that found by 
(2008") in the giant radio galaxy 6C0905-h39 and by 
[2009) for HDF130 in the Chandra Deep Field North. 



The combined pressure produced by relativistic electrons and 
the magnetic field is P/ke = 5.2 x 10^ Kcm~^, which is com- 
parable to that near the centre of a group of galaxies, but much 
larger than that of the IGM (- 0.3 K/cm"^). Therefore whether 
the lobes continue to expand depends on the environment. 



4 MRC SOURCES 

We also obtained observations of giant r adio galaxies MRC 2216- 
206 {z = 1.148, linear size 730 kpc dKaoahi et al lll998h) and 
MRC 0947-249 {z ^ 0.854, linear size 530 kpc (iKapahi et aP 
'1998')) with XMM-Newton on 2008 May 03 and May 10, respec- 
tively, with corresponding exposure times of 22.6 ks (8.8 ks]3 and 
17.0 ks (10 ks). In the 0.3 - 10.0 keV band, there are 361 ± 19 
counts (of which 181 =b 6 are background, with 47 ± 19 sof{fl pho- 
tons) for MRC 2216-206 and 195 ± 14 counts (of which 126 =b 6 



^ Good time remaining after filtering to remove flares and taking dead-time 
intervals into account 
3 0.3 - 1.5 keV 



are background, with 46 =b 11 soft photons) for MRC 0947-249. 
Whereas at higher energies the X-ray images are more consistent 
with being point-like sources near the location of the correspond- 
ing optical host galaxy, indicated by crosses (Fig. |5] and [6]), both 
sources show extended soft X-ray emission below 1.5 keV, which 
is plausibly due to ICCMB. 

We note that the soft emission in MRC 2216-206 appears to 
extend beyond the location of the South-East radio hotspot. This 
could indicate an earlier outburst resembl ing the dou ble-double ra- 
dio galaxy phenomenon ( Schoenmakers et al ][2003) but with the 
outer/older lobes only being detectable via ICCMB. However this 
extra emission, which consists of 30 counts in the 0.3 — 1.5 keV 
band, may also be due to a background AGN (at the corresponding 
flux level of 6 X 10"^^erg s"^ cm"^ in the 0.5 - 2.0 keV band 
for a photon index of F = 2.0 for these 30 counts, we estimate 
a 30 per cent probability of a coincident background source in the 
regions of interest around our three sources using the source counts 
of Hasinger et al. 2001). We did not include the region beyond the 
radio hotspots in our analysis. 

Spectral analysis of MRC 2216-206 was carried out exactly 
as for 3C 469.1, since a single-component power law fit gave an 
unphysically hard photon index of F = 0.5 =b 0.35. The two- 
component model yielded an extra host absorption of 3.8 =b 1.8 x 
10^^ cm~^ and a best fit energy index, F = 1.61q;4, consistent 
with a radio synchrotron index of a = 1 derived from NED archival 
radio photometry. The unabsorbed 2 — 10 keV luminosity of the nu- 
cleus is 1.3 X 10^^ erg s~^. The unfolded spectrum, along with the 
best fit model and residuals is shown in Fig.|7] 

If we assume that the observed X-ray s are inverse-Comp ton 
CMB photons, we can follow equation 4 of lErlund et al.l (l2006h to 
estimate the total energy in relativistic particles from the X-ray flux 
alone: 



L44 



4 7e(l+^)4 



10^^ erg, 



(5) 



where L44 is the X-ray luminosity in 10^^ erg s ^ and 7e is the 
typical Lorentz factor of the electrons up- scattering the CMB to X- 
rays. Setting 7e = 10^ gives 8^ = 5.7 x 10^^ erg for MRC 2216- 
206. This is a lower limit to the total particle energy in the lobes of 
this object. 

There were insufficient counts to fit the spectrum of 
MRC 0947-249. Instead, the soft photon count rate (4.6 x 10"^ 
s~^) was used along with a power-law model and an assumed pho- 
ton energy index of F = 2.0 to estimate the flux of the extended 
emission with PIMMS. The luminosity and total energy in relativis- 
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Source 


RA 


Dec 


z 


countrate 




Unabsorbed 


L44 






(J2000) 


(J2000) 




[cts/s] 


[10^1 cm-2] 


Flux 






3C 469.1 


23h55m23.3s 


-H79d55m20s 


1.336 


1.93 X 10-2 


1.38 


3.58 


2.56 


8.1 


MRC 2216-206 


22hl9m44.2s 


-20d21m31s 


1.148 


2.09 X 10-2 


0.25 


2.93 


1.63 


5.7 


MRC 0947-249 


09h49m52.9s 


-25dllm40s 


0.854 


4.50 X 10-3 


0.46 


1.22 


0.45 


2.9 



Table 2. Observed X-ray luminosities for the sources in this paper. A^h is the Galactic neutral hydrogen column density along the line 
of sight, calculated using PIMMS. Columns 7 and 8 show the observed X-ray flux in 10- erg cm-^ s-^ and the X-ray luminosity in 
10^^ erg s-i (2 — 10 keV, rest frame) respectively, both corrected for Galactic absorption. Column 9 shows the total energy in relativistic 
electrons in 10^^ erg, estimated using equation (|3j for 3C 469.1 and equation (|5) for the MRC sources. For 3C 469.1 and MRC 2216 the 
quoted values of flux, luminosity and energy are for the lobe component as determined from spectrum-fitting, while for MRC 0947 they 
have been calculated from the observed soft (0.3 — 1.5 keV) X-ray count rate of 4.6 x 10-^ s-^ using PIMMS. 
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Figure 7. Unfolded XMM pn spectrum of MRC 2216-206. The model fitted 
is described by equation [T] Note that the y-axis is equivalent to uFi, . The 
solid red line depicts the extended ICCMB component and the green solid 
line the absorbed nuclear component. 

tic particles thus determined are given in Tabled together with the 
corresponding values for 3C 469.1 and MRC 2216-206 determined 
from spectrum-fitting. 



5 CONCLUSIONS 

We have detected extended X-ray emission in the giant FRII radio 
galaxy, 3C 469.1, which extends up to the northern radio hotspot. 
The spectrum is well fit (x^ = 13.9 for 17 d.o.f.) by a power 
law with a photon index consistent with the radio spectral index. 
Energy equipartition arguments agree with the magnetic field es- 
timates made for 3C 469.1 based on the X-ray flux inferred to be 
ICCMB, provided 7min is in excess of 250. 

Further, we have presented XMM images of two more radio 
galaxies at cosmological distances, which show extended X-ray 
emission attributable to inverse-Compton CMB. The nuclei of all 
these sources have quasar-like luminosities and appear to be heav- 
ily absorbed powerful AGN, while the radio-emitting lobes have 
relativistic particles with similar total energies, ranging between 3 
and 6 x 10^^ erg. These values are enhanced if there is significant 
energy in protons. 

It appears that ICCMB X-ray emission is indeed observable 
from giant radio galaxies at redshift z ^ 1. However, the cospatial 
nature of the X-rays with radio observations remains an important 
caveat in the determination of the magnetic field. Further observa- 
tions of these sources with the finer PSF of Chandra would help re- 



solve this issue and observations with greater signal-to-noise would 
provide better constraints on the X-ray spectral index, which is es- 
sential for accurate determinations of 7min and the establishment 
of the nature of the low-energy turn-over in aged AGN lobes. 
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